New M-dwarf debris disk candidates in NGC 2547 
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ABSTRACT 

With only six known examples, M-dwarf debris disks are rare, even though M dwarfs constitute 
the majority of stars in the Galaxy. After finding a new M dwarf debris disk in a shallow mid- 
infrared observation of NGC 2547, we present a considerably deeper Spitzer-MIPS image of the 
region, with a maximum exposure time of 15 minutes per pixel. Among sources selected from 
a previously published membership list, we identify nine new M dwarfs with excess emission at 
24 fim tracing warm material close to the snow line of these stars, at orbital radii of less than 
1 AU. We argue that these are likely debris disks, suggesting that planet formation is under way 
in these systems. Interestingly, the estimated excess fraction of M stars appears to be higher 
than that of G and K stars in our sample. 



Subject headings: circumstellar matter — infrared: stars 
(NGC 2547) — planetary systems: formation 



1. Introduction 

Even though about 80 % of t he stars in the 
Galaxy are M dwarfs ( Lada! 120061 ). only very few 
of these are known to harbor debris disks. This 
is even more surprising since many M dwarfs 
do show signs of protoplanetary disks in ear- 
lier stages of evolution wi th similar frequency 
as more massive stars (e.g 



Andrews fc Williams 



20051 iMuzerolle et al.ll200a lLada et al.ll2006l ) and 



we now know seven M dw arfs that harbor extraso- 
lar planets (most recentlv rEndl et al]|2008i see ref- 
erences therein). While protoplanetary disks con- 
tain large amounts of primordial gas and dust out 
of which planets form, debris disks are older and 
entirely made of collisionally evolving planetesi- 
mals and dust particles. Debris disks, thus, are 
secondary products of the planet formation pro- 
cess and their existence in a given system implies 
previous formation of protoplanets. Compared to 
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open clusters and associations: individual 



optically thick protoplanetary disks, debris disks 
are fainter and optically thin. For a recent ob- 
servational and th eoreti cal re view of debris disks, 
(|2007l ) and iMoro-Martin et al. 



sec 



Meyer et al. 



(120081) . 

Since the cold dust in debris disks is best stud- 
ied at far-infrared (FIR) or submillimeter wave- 
lengths, the Spitzer Space Telescope, given its 
unprecedented FIR sensitivity, has revolutionized 
the research on the role of debris disks in star 
and planet formation. Most of the Spitzer stud- 
ies of debris disks have been carried out using 
the Mult iband Imag i ng P hotometer for Spitzer 



(MIPS, iRieke et all 120041) . operating at wave- 



lengths of 24 /im, 70 /an, and 160 fim. 

Prior to our study, only six debris disks around 
M dwarfs were known. Early studies based on 
data obtained by IRAS, the Infrared Astronomi- 
cal Satellite, reported far-infrared excess emission 
towa r ds a small number o f field M stars (ITsikoudi 
1988[ iMullan et al] Il989t iMathioudakis fc Dovle 
1991), but few of these results could be confirmed 



subse quently (e.g., ISong et all l2002t iRiaz et al 



2006). The first firm detection of an excess was 



the disk of AU Mic (JKalas et al.ll2004 |Liu| 12004), 



a member of the f3 Pic moving group with an 
age of ~12 Myr. The discovery observation was 
made with a stellar coronograph detecting scat- 
tered light in near-infrared broad-band filters. 
AU Mic was also found to ha ve excess em ission 
at 70 fim, but not at 24 /mi (jChen et all 12005) . 
Subsequently, in subm il limete r radio observations 
at 850 fim, iLiu et all (|2004l ) identified GJ 182, 
a member of the Local Association Group with 
an age of «50 Myr, as harboring a debris disk. 
GJ 182 does not show a far-infrared excess in 



the Spitzer-MIPS study of IChen et all (2005). In 
Spitzer-MIPS data of NGC 2547 a clu ster at an 
age of 30-40 Myr, lYoung et all (|2004 ) detected 
the third M star debris disk (source 23 in our Ta- 
ble [2l see also Teixeira et al, in prep.). TWA 7 
was the fourth object reported; it is a member of 
the TW Hya association with an age of < 6 Myr 
and excess emission reported at 24 /i m, 70 /mi, 
and s ubmillimeter radio wa velengths ( Low et all 
120051 iMatthews eFaTlteOOTl) . One field M dwarf 
(GJ 842.2) was found to have a debris disk (i n 
submilli meter observat i ons, Lestrade et all 120061 ) . 
Finally, iGorlova et all (|2007l ) identified one addi- 
tional M dwarf debris disk in NGC 2547 at 24 /mi, 
listed as source 22 in our Tabled] 

Only three sources (TWA 7, source 22, and 
source 23) have been identified as far-infrared ex- 
cess sources at 24 /tm, while the other sources 
have been found at either 70 /im (AU Mic) or 
submillimeter radio wavelengths. Since material 
of quite different distances from the central star 
is probed with these two techniques, debris disks 
identified in the mid-infrared and in submillime- 
ter radio are not necessarily directly comparable 
to one another. 

Studying a sample of 62 field M dwarfs es 



timated to be older than 1 Gyr, Gautier et al. 



(J2007Q did not find any indication of excess emis- 
sion at 2 4 /tm or 70 /xm indica t ive of debris disk s 



see also IPlavchan et all 120051 iRiaz et al. 2006) . 



In con trast to these Spitzer results, Hestrade et al 



( 20061 ) discuss two sets of (sub-)millimeter data, 
including previously published observations, and 
find three detections of excess emission in a sam- 
ple of 23 M dwarfs with ages ranging from 20- 
200 Myr, corresponding to an excess fraction of 
13_ 8 %, which may not be significantly different 
from the disk fractions of earlier spectral types. 
After identifying an M dwarf debris disk in ear- 



lier Spitzer observations of NGC 2547 (source 23, 
see above), we performed ten times longer Spitzer- 
MIPS 24 /im observations of the same region in 
order to find out whether this was a singular ob- 
ject or part of a larger population. In the remain- 
der of the introduction, we summarize the current 
knowledge about NGC 2547 before describing the 
observations in Section [3] We present the results 
in Section |H including remarks on the 70 /im and 
160 /tm data as well as variability at 24 /tm, dis- 
cuss these results in Section [5] and close with a 
summary in Section [BJ 

2. Background: Cluster Parameters 



NGC 2547 i s a young open cluster at a distance 
of 4331^ pc (jRobichon et al.lll999l ) and an av- 
erage visual extinction of only Ay — 0.18 mag 
( Clarial ll982J). It has an ap proximate age of ~ 
30-40 Myr (JLvra et alj200a who fit four different 
evolutionary tracks to UBVRI photometry, three 
of the m yielding the older age). iNavlor &: Jeffries! 
([20061 ) determine an age of 38.5+^ Myr (and 
a lower distance of 36lj^g 9 pc) from a fit to a 
color-magnitude diagram, improving on an ear- 
l ier age estimate i n the range of 20 - 35 Myr 



(Naylor et al. 



,. 2002, based on BVI photometry). 
Young et al.l ( 20041 ) discuss the disk population 
based on Spitzer imaging data, including one 
epoch o f the MIPS da t a pre sented in this paper, 



as does IGorlova et all (|2007l ). Both of them use 



the data we mark as "epoch 1" . 



Gorlova et all (|2007l ) provide a membership 



list based on a dee p RIZ survey carried out by 
Jeffries et all (|2004l ). While the latter paper dis- 



tinguishes empirical candidate me mbers derived 



using different evolutionary tracks, IGorlova et al. 



(|2007l ) use photometric, spectroscopic, and proper 
motion information for different subsamples to 
produce a single list of 806 "highly probable mem- 
bers" . However, 64 of these sources appear not to 
warrant such a classification because they are lo - 



cated in the two fields that Ijeffries et all |200J) 



used to check for galactic contamination, far from 
the actual cluster (outside the coverage of Fig. [1]) . 
Therefore, the total number of "highly probable 
members" is 742. Twelve sources in the full list 
do not have 2MASS-K mag nitudes, althou g h they 

(2007) 



are detected in J band. Gorlova et al 



also provide near- infrared spectra for 89 sources. 
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Fig. 1. — MIPS-24 jUin map of NGC 2547. The black line delineates the "core area" with at least two 
thirds of the maximu m exposure while the white line contains the "highly probable members" as defined by 
Gorlova et al.l ([2007J). The image scale is logarithmic and the pixel size is 2.45". 



Irwin et alJ ( 20081 ) determine rotation periods at 
optical wavelengths for 176 photometrically se- 
lected candidate members. 128 of these have 
counterp arts in the "h i ghly probable members" 
listed bv lGorlovaet alJ (|2007l ). 



3. Observations and Data Reduction 

We present a deep Spitzer-MlPS observation 
of the central part of NGC 2547. The combined 
dataset that we are discussing here consists of all 
MIPS maps taken of NGC 2 547 to date; it i nclude s 
the one map analyzed by iGorlova et al.l ( 2007T ). 
We took ten additional maps, but due to an in- 
strument failure during observation, one map was 



lost. Thus, we discuss a total often different maps 
that can be grouped into three epochs (Table [1]). 
All but epoch 1 also have 70 ^m data. Additional 
160 //m data was taken in the three observations 
of epoch 2. 

The entire 24 /jm dataset consists of 21127 
"basic-calibrated data" (BCD) files, processed 
with pipeline S14.4.0. These data were reduced 
using the Spitzer M Osaicker and Point sourc e 
Extractor (MOPEX. iMakovoz fc Marleaul l2005h . 
version 030106. In order to produce a single deep 
map out of single maps with different background 
levels, the MOPEX task overlap.pl was used on 
subsets of the full dataset after flatfield correction. 



The subsequent data reduction steps including 
median filtering, outlier detection and masking 
before coadding and combining were carried out 
on these subsets of BCDs. The resulting map is 
shown in Fig. Q] The map shows that the region 
is largely devoid of extended emission, as faint ex- 
tended features can only be seen in the logarithmi- 
cally stretched image. Subsequently, the Spitzer 
Astronomical Point Source Extraction (APEX) 
software was used for photometry on the resulting 
maps, applying point-response-function (PRF) fit- 
ting with local noise estimation and a cutoff based 
on the signal-to- noise ratio (SNR). Prior to the ac- 
tual photometry, a sample PRF was determined 
from several well-detected stars in the field. A 
number of objects were detected in several of the 
overlap-corrected map chunks which themselves 
are overlapping. In these cases, photometry was 
taken from the map with the higher exposure map 
value in order to ensure optimal coverage. The 
photometric accuracy was checked by subtracting 
the fitted sources from the original maps in order 
to create residual maps. Only two sources have 
residuals that are clearly off (even negative); both 
are bright sources within extended emission (see 
also discussion below). Corresponding steps were 
taken for the analysis of the 70 /jm and 160 /zm 
data. For comparison, point-response function 
photometry was also performed on the newly re- 
duced archival IRAC data. 



Table 1: NGC 2547 MIPS 24 /im observations 



Epoch 


UT start date & time 


further obs. 


1 


2004-01-28 19:17:57.6 


- 


2a 


2006-05-08 23:01:46.6 


70, 160 /im 


21) 


2006-05-09 07:29:42.3 


70,160/im 


2c 


2006-05-09 12:02:06.3 


70, 160 /im 


3a 


2006-12-01 22:01:58.8 


70 /im 


3b 


2006-12-02 01:03:24.0 


70 (im 


3c 


2006-12-03 21:36:33.4 


70 /im 


3d 


2006-12-04 00:38:00.7 


70 /im 


3e 


2006-12-04 03:44:48.9 


70 /im 


3f 


2006-12-04 06:46:16.2 


70 /im 



4. Results 

4.1. 24 /im data 

4-. 1.1. Source detection 

At a 3<t detection level, a total of 27042 sources 
is detected. This is our full sample with all detec- 
tions in the full coverage area, including outlying 
areas which were covered only once. The effective 
exposure time per pixel in this full dataset varies 
by a factor of ten across the map. In order to have 
a more homogeneous dataset, we only use detec- 
tions with at least 2/3 of the maximum effective 
exposure time per pixel (t max = 901.6 s) in a sec- 
ond sample. This second sample, which we will 
refer to as the core area, focuses on the central re- 
gion of NGC 2547; it includes 6390 sources. The 
distribution of 24 /im magnitudes in this core area 
peaks at a value of 11.6 mag which compares to 
a peak at 10.8 mag in the epoch 1 data used by 
Gorlova et al.1 (|2007l ). Every single map roughly 
has 10% of the total maximum effective exposure 
time per pixel. 

As a first step, the 24 /im detections were corre- 
lated with the 2MASS point source catalogue since 
all candidate members should be NIR-detected. 
There are 1945 sources from the full list that have 
2MASS counterparts within l'/5. This radius cor- 
responds to 3 times the half-width half-maximum 
of a Gaussian distribution fitted to the angular 
separations between the 2MASS and the MIPS 
data; it was used for all subsequent correlations. 
Fig. [2] shows that there is a good correlation be- 
tween 24 /im and K band magnitudes down to 
Kk, 12. The scatter of 24 /im magnitudes for 
sources with a given K magnitude of K > 12 mag 
indicates that these are likely galaxies instead of 
single stars since 2MASS photometric errors in K 
increase significantly only beyond K = 14 mag. 

In fact, only 646 sources out of the 6390 sources 
in the core area have 2MASS counterparts, in- 
dicating that the overwhelming majority of the 
sources we detect at 24 /im are galaxies that are 
too faint in the K band to be detected by 2MASS. 
Out of the 6390 sources, 208 have counterparts 
in the Gorlova membership table. Correlating all 
27042 sources in our entire dataset to the Gorlova 
list yields 265 counterparts. This means that the 
white polygon in Fig. [T] encompasses 265 detec- 
tions while the black polygon contains 208 detec- 
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Fig. 2.— MIPS-24 /am vs. K band (2MASS) pho- 
tometry for the 1945 MIPS sources with 2MASS 
counterparts; error bars are nominal la. A few 
sources have large error bars and some of the 
faintest 2MASS sources shown do not have pho- 
tometric errors. 
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Fig. 3. — Comparison of the Rc band magnitudes 
from ljeffries et al.l ( 20041 ) with R band magnitudes 
from USNO-B1 for the 208 highly probable mem- 
bers detected at 24 /im in the core area. For 
sources brighter than 12.3 mag (dashed line), sat- 
uration effects become important, a nd we use the 
USNO R photometry instead of the Ijeffries et al. 
(J2004I ) photometry. 
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Fig. 4. — Estimate of completeness in the core 
area. The solid-line histogram shows the R-K 
distribution of the 554 candidate members lo- 
cated in our core area (and having R-K pho- 
tometry), the dashed histogram shows our 208 
detections at 24 /im, and the dotted histogram 
shows the 25 detect ions at 24 ^m reported by 
Gorlova et al.l IJ2007I ). The indicated spectral 
types from B stars to about M6 are based on 
R-K and are for luminosity clas s V, u s ing in- 
trinsic colo rs fromlBessell fc Brettl (119881) . iBessell 
(J1991I ). an dlKenvon fc Hartmannl (|l995l ) with up- 
dates from lWinklerl (J1997I) . 



tions. iGorlova et al.l (J2007I ) extracted 74 candi 



date members with 24 /zm counterparts (from the 
map we list as Epoch 1); 64 of these are detected 
in the core area although it actually covers 71 of 
them. The seven missing sources either did not 
pass our SNR selection criterion or the angular 
distance between a detection and the catalog po- 
sition was too large. Th e median magn i tude d iffcr- 
ence at 24 /zm between IGorlova et al.l (|2007l) and 
our dataset is 0.04 mag, including five cases where 
the difference is larger than 0.5 mag (these are 
the bright sources within extended emission; the 
maximum difference is 1.5 mag). 

The following discussion focuses on the core 
area. For the near-infrared photometry, we make 
use of the 2MASS and IRAC photometry from 
Gorlova et al.l (J2007I ). The deepest optical data 
come from t he R r I r Z survey carried out by 
Jeffries et al.l ( 20041 ). We therefore focus on the 
R-K color baseline. 

Fig. [3] shows a direct comparison of photo- 
graphic R band magnitudes from the USNO-B1 



catalog with those from I Jeffries et al.l (|2004j ). Be- 
cause of R c band saturation problems for the 
brightest sources in this sample (i? c < 12.3 mag), 
we use R magnitudes from the USNO-B1 cata- 
log (JMonet et al.l 120031 ) instead for those bright 
sources. The median differ e nce in R between 
USNO-B1 and Ijeffries et all (J2004 within the 



of photometric members are needed since there is 
a lot of contamination, espec ially due to the giant 



208 candidate members is 0.13 mag, but it is 
considerably higher for the brightest sources (up 
to 3.9 mag). The photometric accuracy of the 
USNO-B1 catalog is given as 0.3 mag. By combin- 
ing these two sets of photometry measurements, 
we can make use of reliable R band photometry 
for all 208 sources in our core sample, at the price 
of varying accuracy. For comparison, th e R band 
photometry used bv lGorlova et al.l (|2007l ) does not 
include information for 47 out of the 208 sources. 

There are a few exceptions to our R < 12.3 mag 
selection rule. Due to a binary source fitted as a 
single source in US NQ-B1, two sources with R < 
12.3 mag only have I Jeffries et al.l (j2004h photom- 
etry which we therefore use as best choice. There 
are nine other sources in the sample which have 
unreliable USNO-B1 photometry for the same rea- 
son. Also, two obviously faulty R band magni- 
tudes in the core sample (R > 28 mag) were re- 
placed by the USNO values. In total, we make 
use of USNO-B1 photometry for 53 out of the 208 
sources. 

In Fig. [U we show the distributions in R-K 
of the 562 "highly probable membe rs" with R-K 



photo metry in our core area (using iJeffries et al. 



120041 and USNO-B1 photometry), our 208 detec 
tions at 24 //m, as well as those det ections within 
our co re area that were reported bv lGorlova et al. 
(|2007l ). These authors report R-K for only 26 out 



of their 74 detections at 24 /im. For our analysis, 
we use the same combination of USNO-B1 and 



Jeffries et al.l (2004) photometry for these sources 



that we also use for our own data to provide R 
mag nitudes: seven t y-one of the 24 fim sources 



from lGorlova et alj (|2007l ) are located in our core 
area and seventy now have reliable R band mag- 
nitudes. While by far not detecting all candidate 
mem bers, we detect a con siderably larger fraction 
than lGorlova et al.l ( 20071 ). Fig. [^indicates a com- 
pleteness dropping to ^50% in the K stars. How- 
ever, although not complete, we still get a good 
sampling of K and M stars. It becomes immedi- 
ately clear that more spectroscopic confirmations 



Jeffries et all (|2004f ) 



branch at R-Kfa2.5. In fact, 
identify the region with 0.5<-R-I<0.7 as particu- 
larly prone to contamination. There are 167 out of 
562 candidate members in our core region that fall 
into this range which corresponds to a full range 
in R-K of 1.84 to 2.68. However, as will become 
clear later, this does not affect the discussion of 
the M stars. 

Looking at the entire dataset, only 120 out of 
265 highly probable members detected at 24 /im 
have available V-K colors, ranging from -0.36 to 
5.73. In the core area, this ratio is only 111 out 
of 208. Since, with the above procedure, we have 
R-K for all 208 sources of interest, we used R-K 
colors instead of V-K as a proxy for spectral type. 

4-1.2. Excess sources 

We used two different criteria to identify 
sources w ith excess emi s sion a t 24 /xm. Initially, 
following iGorlova et al.1 (|2007l ). we used the 3ct 
locus of the Pleiades stars in a K-24 vs. R-K dia- 
gr am of our data . The Pleiad es locus was derived 
bv lGorlova et al.l (|2006ll2007t ) . However, since the 
membership list for NGC 2547 is not yet defini- 
tive, we also looked at the distibution of all sources 
detected at 24 /mi in the same area for which R 
and K photometry is avail able, including 1 4 "les s 



Gorlova et all (|2007h . 



plausible" members from 
To crudely exclude extragalactic sources, we only 
use candidate non-members with K < 12 mag (see 
above), yielding 275 sources in the core area. We 
characterize their location in the color-color dia- 
gram by determining the standard deviation in 25 
R-K bins of equal sample size, clipping the most 
extreme point in every bin. The resulting color- 
color exclusion zone is defined as the 3tr variation 
in every bin. Note that for the i?-band photom- 
etr y, we use the same combination of USNO-B1 
and lJeffries et al.l ( 20041 ) photometry that was also 
used for the members (see above). 

According to these criteria, 54 sources have sig- 
nificant excess emission compared to the Pleiades, 
i.e., their 3er photometric error bars lie completely 
outside the 3er locus of the Pleiades. Twenty- 
six of these sources also are outside the color- 
color exclusion zone (including sources embedded 
in bright extended emission) . In the following, we 
use the latter as the more conservative criterion 
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Fig. 5. — Plots of K-24 /jm vs. R — K for a) 275 24 /xm-detected sources in the cor e area and not in the 
"highly probable member" list (14 matches with the "less plausible" members from iGorlova et alj ( 2007 ) 
are marked by 'o', but only a few are slight outliers), b) 208 24 /im-detected highly probable members in 
the core area, c) the 54 significant excess sources according to the Pleiades criterion, and d) the 24 excess 
sources outside the color-color exclu sion zone; boxes den ote the four significant excess sources (using our 
exclusion zone criterion) identified bv lGorlova et al.l ( 2007) in the cor e area. The solid lines indicate the locus 
of the Pleiades and respective 3a limits, as used by IGorlova et al.l (|2007r ). and the dashed lines delineate 
the 3a color-color exclusion zone (see text). A measure of completeness, the dotted line shows the mean 
of the distribution o f 554 "highly probable members" sources with R and K magnitudes in the core area 
(jGorlova et al.l 120071 ). assuming a limiting 24 /xm magnitude of 11.6 mag. Sources with large error bars 
in R — K have R band photometry from USNO-1B (see text). The 'x' symbol in panels b), c), and d) 
m arks five source s that are located in bright extended emission and were already designated as "halo-like" 
by IGorlova et al.l ( 2007 ). The error bars show the propagated la errors. The spectral types indicated in 
panel d), from B stars to about M6, are based on R-K and are for luminosity class V. 



to define excess sources. The results of the search 
for excess sources are summarized in Fig. [5] and 
in Table [2] Panel a) of Fig. [5] shows the defi- 
nition of the color-color exclusion zone by plot- 
ting all 275 candidate non-member sources within 
the core area having 24 /im detections as well 
as R and K magnitudes (with K < 12). Panel 
b) shows all 24 /im-detected "highly probable 
members" (jGorlova et al.l 120071 ) within our core 
area, with the near-infrared photometry explained 
above. The lower panels show the significant ex- 
cess sources when using the Pleiades (panel c) 
and the color-color exclusion zone (panel d) cri- 
teria. Both Pleiades locus and exclusion zone are 
shown in all plots. Panel d) additionally shows 
the previously identified four excess sources within 
the core area (using the exclusion zone criterion) 
(IGorlova et al. 20071) as circl es. A fifth excess 
source in IGorlova et al.l <|2007l ). at R-K= 1.95 mag 
and K-24=1.91 mag, is not detected in our data 
(and not shown in the plot) because the positional 
difference between the 24 /im source and the NIR 
counterpart is too large (3.6" in our case). Since 
the 24 /zm magnitudes are nearly identical, this is 
probably the same source, however. 

We estimate sensitivity across R-K by assum- 
ing a limiting magnitude at 24 /im and calculat- 
ing the distribution in K-24 fun vs. R-K of the 
554 "highly probable members" located in the core 
area that have suitable photometry (eight sources 
lack 2MASS-X photometry). We then determine 
the standard deviation in K-24 /im for 46 bins con- 
taining nearly the same sample sizes, with a fixed 
value at 24 /im of 11.6 mag (see above), thus char- 
acterizing the distribution per R-K bin. Also this 
analysis takes into account USNO-B1 photometry 
for the brightest sources. Panel d) of Fig. [5] shows 
the mean of this distribution. It becomes clear 
that for late M dwarfs, we can only detect ex- 
cess sources, whereas for earlier spectral types, we 
can also detect photospheres, although the detec- 
tions of K and M dwarfs are not complete. In fact, 
most of the M dwarfs that we detect at 24 /im (14 
out of 25) do not show significant excess emission. 
To corroborate this result, we show the SED of a 
spectroscopically confirmed MO member that we 
detect as a source with a significance of 4.6 a at 
24 /im at presumably photospheric emission levels 
in Fig. [6J 



24 /im detections as well as R and K photome- 
try, of which nine have excess emission compared 
to the Pleiades and six have excess emission with 
respect to our color-color exclusion zone. Only one 
of these sources is not in our core area. Another 
source from within these six sources is 3.6" from 
the nearest 2MASS source in our data, so we do 
not list it. The remaining four previously identi- 
fied excess sources are marked as circles in panel 
d) of Fig. El 

As remarked by IGorlova et al.l ( 2007 ). five of 
the 24 /im-brightest sources covered in both their 
data and our core area are embedded in bright 
extended emission. Four out of these five techni- 
cally have excess emission, even according to the 
color-color exclusion zone criterion, but the ex- 
tended structures are too large to represent debris 
disks. We compared the the PRF-fitting photome- 
try with aperture photometry and found that the 
the former appears to reasonably reflect the ac- 
tual source brightness in the central areas while 
the curve of growth of the latter does not con- 
verge. The photometric results for these bright 
sources with extended emission have to be inter- 
preted with caution. 

Even though we are dealing with relatively low 
numbers, it is interesting to look at the excess frac- 
tion per spectral type, as derived from R-K colors, 
which we plot in Fig. [7] While the contamination 
of the membership list makes it difficult to state 
exact fractions, this problem is mostly constrained 
to the K stars (see Fig. [4j, where many of the 
"highly probable members" may be background 
giants. Even though our data are not complete for 
M dwarfs, 11 of 26 significant excess sources are 
M dwarfs, corresponding to a lower limit of 4.9% 
of all "highly probable" member M stars in the 
area, or 44% of the M dwarfs we detect at 24 /im. 
Four of these 11 M dwarfs are spectroscopically 
confirmed members (Table [2]) . The excess rate is 
is of particular interest given the minimum in the 
G stars, where we detect just a single significant 
excess source. 

Warmer protoplanetary disks closer to the cen- 
tral star would be recognizable by their mid- 
infrared excess emission at shorter wavelengths, as 
probed by the Spitzer-IRAC bands. Referring oth- 
erwise to the discussion of the available IRAC data 



Gorlova et al.l (J2007I ) report 26 sources with 



Young et all (|2004h and IGorlova et alJ (120071 ). 



who found an excess fraction of <1% in the IRAC 



bands, we note in particular that there is virtu- 
ally no sign of excess emission in the IRAC bands 
towards the M dwarfs with 24 /im excess emis- 
sion (eight out of eleven are covered by the IRAC 
data; see also the discussion of the SEDs). The 
only excess de t ection is source 23, listed as ID 7 in 
iGorlova et al.l (2007J), which has an 8<r excess at 
8/iin. 

4.2. 70 /im and 160 /im data 

After data reduction, we performed PRF-fitting 
photometry on the 70 /im and 160 /im data using 
APEX at a 5a detection limit. In an image built 
from the filtered 70 /im BCD data, 390 sources are 
detected. Within search radii of 2" and 4", only 
one of the "highly probable members" is clearly 
detected, with an SNR of 15; this is source 15, also 
showing the strongest excess emiss ion at 24 /im, 



i ts 70 fim magnitude is 2.76 mag. IGorlova et al. 
( 20071 ) concluded that it may be an equal-mass 
binary late K dwarf (their source ID 9). For all 
other sources, upper limits from aperture photom- 
etry were obtained, but they do not constrain the 
spectral energy distributions of these sources and 
thus are not discussed further. 

At 160 /tm, only the extended emission features 
are recognizable, otherwise none of the "highly 
probable members" are detected. Also here, up- 
per limits derived from aperture photometry do 
not constrain the source SEDs. 

4.3. Variability 

The fact that the 24 /im observations were car- 
ried out in three distinct epochs offers the possi- 
bility to study variability on timescales of months. 
Looking for large variability, PRF-fitting photom- 
etry was run separately on all ten maps. Then, by 
cross-correlation, for every source the number of 
detections in these ten maps was determined. As a 
(rather crude) measure of variability, the standard 
deviation of the set of photometry measurements 
was calculated for every source. 

At a 3a detection level, 38 highly probable 
members are detected in all ten 24 /im maps. As 
a measure of variability, we determine for every 
source the ratio of the maximum (presumed) vari- 
ability amplitude divided by the noise level. When 
merging the photometry into one value for each 
epoch, only four out of 38 sources have a variabil- 




Fig. 6.— SED of 2MASS J08095644-4922109, a 
co nfirmed M0 membe r of NGC 2547 (source 37 



Jeffries et al.l 120031 ) without significant 24 /im 
excess emission. The fit uses the Phoenix 



model (JHauschildt et al.lll999aH bl). with log(g)=5 
and solar metallicity; the derived temperature is 
T=4100 K. 
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Fig. 7. — MIPS-24 /im excess sources as a function 
of spectral type. The uppermost histogram shows 
all 554 "highly probable members" with R-K in- 
side our core area, the next histogram shows the 
208 sources we detect at 24 /im, while the lowest 
histogram shows the 26 significant excess sources 
outside of the color-color exclusion zone. The ex- 
cess fractions are given as a percentage of the 554 
candidate members in the area. No percentage is 
given for the B stars because some of them arc 
excess sources with extended emission that clearly 
are not debris disks (see text). The excess fraction 
for K and M stars are marked with an asterisk be- 
cause the data probably are not complete for these 
spectral types (see text). 



Table 2 
mips-24 flm excess sources outside the color-color exclusion zone, sorted in r-k 



2MASS 



K 

[mag] 



R-K 

[mag] 



[mag] 



K-24 
[mag] 



err. 
[mag 



24^tm 
[mag] 



err. 
[mag 



comments 



SNR 
24/jm 



1 


J08102058-4914144 


6.83 


-0.33 


0.33 


2.82 


0.03 


4.01 


0.00 




B3III/IV, X, G, extd 


17.1 


2 


J08101607-4902058 


7.65 


-0.29 


0.32 


1.06 


0.02 


6.59 


0.00 




B5II, X, G, extd 


54.0 


3 


J08102725-4909509 


8.06 


-0.14 


0.33 


3.45 


0.03 


4.61 


0.00 




B6V, X, G, extd 


119.9 


4 


J08095753-4908202 


8.95 


-0.01 


0.32 


1.89 


0.02 


7.06 


0.00 


2 


AOVn, G, extd 


89.2 


5 


J08100607-4914180 


9.66 


0.17 


0.32 


1.53 


0.03 


8.13 


0.01 


2 


A1V 


17.8 


6 


J08092668-4914371 


10.14 


0.32 


0.32 


1.35 


0.03 


8.79 


0.01 


2 




46.2 


7 


J08100841-4900434 


9.55 


0.39 


0.32 


1.19 


0.03 


8.37 


0.01 




A3V 


77.0 


8 


J08111134-4904442 


10.18 


0.59 


0.32 


0.73 


0.04 


9.46 


0.02 






12.4 


9 


J08084571-4923473 


9.95 


0.63 


0.32 


1.06 


0.03 


8.89 


0.01 






42.5 


10 


J08110323-4900374 


10.15 


0.72 


0.32 


1.20 


0.03 


8.96 


0.01 






53.0 


11 


J08093815-4918403 


9.67 


0.92 


0.33 


1.00 


0.04 


8.67 


0.01 


2 


A9V, G ID 2, X 


39.2 


12 


.108090250-4858172 


11.36 


1.37 


0.02 


3.92 


0.02 


7.44 


0.00 


2 


G ID 8 


162.1 


13 


J08101799-4923502 


11.65 


1.75 


0.03 


0.72 


0.09 


10.93 


0.07 


2 


X 


9.2 


14 


J08102955-4922094 


12.35 


2.22 


0.03 


1.17 


0.11 


11.18 


0.09 






6.0 


15 


J08101691-4856291 


11.12 


2.59 


0.03 


6.87 


0.02 


4.25 


0.00 


2 


G ID 9, 70 fim 


1500 


16 


J08095330-4918132 


12.85 


3.35 


0.04 


1.37 


0.15 


11.48 


0.12 


4 


M2, RV, X 


3.9 


17 


J08110007-4904413 


13.05 


3.35 


0.04 


1.37 


0.17 


11.68 


0.14 






3.3 


18 


J08101374-4917112 


13.16 


3.59 


0.04 


1.34 


0.19 


11.82 


0.16 




X. P 


3.1 


19 


J08105210-4921138 


12.52 


3.71 


0.04 


1.18 


0.13 


11.34 


0.10 




X 


6.5 


20 


J08092450-4928575 


12.68 


3.80 


0.04 


1.35 


0.13 


11.33 


0.10 




P 


5.5 


21 


,708094510-4846520 


12.82 


3.94 


0.04 


1.39 


0.15 


11.43 


0.12 






5.9 


22 


J08091770-4908344 


13.85 


4.41 


0.06 


3.12 


0.11 


10.73 


0.06 


1,3 


M4.5, RV, P, G 


6.8 


23 


J08093547-4913033 


14.07 


4.41 


0.06 


3.48 


0.09 


10.59 


0.05 


1,2,3 


M4.5, RV, G ID 7° 


8.5 


24 


J08085407-4921046 


14.76 


4.52 


0.11 


2.77 


0.29 


11.99 


0.20 


1 




4.8 


25 


J08104437-4939001 


13.82 


4.75 


0.06 


4.18 


0.07 


9.64 


0.02 


3 


M5, RV 


25.2 


26 


J08101391-4939481 


14.11 


5.01 


0.09 


3.19 


0.13 


10.92 


0.07 


l d 




7.7 




ir-infrared spectrum in 












(120071), 3) 










a Nc. 


1) lOliveira et al 


.M2003), 


2)IGorlova et al.l 


UefTries 


& Oliveira! J2005I), 4)IJeffries et al 


J 12003). 



b sp ectral type from SIMBAD, 'X' = XMM detect ion within < 25", fr om lJeffries et"al] < |2006D ,'P' = optical peri od from [I 
(2008), 'G 1 — previously identified 24 fira source (Gorlova ct al. 2007), 'RV — confirmed rad. vcl. 



Jeffries & Oliveira 2005), 'extd' — in bright extended emission, possibly unreliable 24 fim photometry 
c discovered bv lYoung et al.l (2004) 
d Lithium absorption with EW L i = 0.56 ± 0.27 A 



. rwin et al.| 

member (Jeffries et al. 2003; 
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Fig, 8. — SEDs fitt ed with disk models from 
Robitaille et al.l ( 20071 . solid line) and the Phoenix 
atmosphere model (dashed line) for the five ex- 
cess sources with R-K> 4. Two sources are lo- 
cated outside the footprint of the IRAC mosaic, 
for them, only the atmosphere model is shown. 
Source numbers are indicated in the upper right 
corner of every panel. 



ity exceeding the respective 5a level, and only two 
sources exceed la. These are 2MASS J08101546- 
4905487, one of the brightest sources in the field, 
varying by 19%, and 2MASS J08085649-4923128, 
varying by 24%; both sources are not among the 
excess sources discussed above. 

5. Discussion 

There are three different possibilities for the in- 
terpretation of the obser ved excess emission (e.g., 
Strubbe fc Chiandl2006[) . It could be due to rem- 
nant primordial material, due to a debris disk with 
continuously replenished dust, or it could be debris 
produced by a recent catastrophic collision of plan- 
etesimals, i.e., debris not in (quasi-) equilibrium. 
The dust in debris disks is produced and replen- 
ished by collisions, and then removed by processes 
such as the Poynting-Robertson (PR) drag, radi- 
ation pressure and stellar winds. Direct radiation 
pressure does not play a role for the low-lumi nosity 
stars that we consider ( Plavchan et al.ll2005l ). The 
balance of dust production and these removal pro- 
cesses determines the amount of observable dust. 

In order to discuss the physical processes, we 
first need to determine which region around a cen- 
tral star is probed by the 24 /im (excess) emission. 
As an estimate, we consider for a given stellar lu- 
minosity the orbital radius of a blackbody which 
is at a temperature that puts the peak of its SED 
at 24 /jm. For a dwarf star with a luminosity 
of i=0.01 Lq, the corresponding radius is only 
0.5 AU. 

With the orbital radius of the dust particles 
constrained, we can now determine the lifetime 
of dust particles under the influence of the PR 
drag, which should yield an upper limit of the dus t 
removal timescale (e.g... iDominik fc Decinl 12003 ). 
The PR lifetime of a dust particle in orbit around 
a central star scales with particle size and den- 
sity as well as the square of its orbital radius and 
is inversely proportional to the luminosity of the 



central star (Burns et al 



19791 ). For micro n-sized 
particles with a density o f p = 2.5 gem -3 (jDivinei 
1993t IChen fc Jurall2001l ) located at a distance of 
1 AU from a central star with a luminosity of 
L = 0.011/0, the lifetime is still on the order of 
1.7 x 10 5 years. This value is ten times longer 
for 10 /zm-sized particles. Since the M dwarfs dis- 
cussed here have an age of 30-40 Myr, the particles 
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would need to be replenished and thus cannot be 
primordial in order to still be observable at that 
age, indicating the debris disk nature of these sys- 
tems and ruling out the first possibility mentioned 
in the beginning of this section. The lifetimes do 
not approach the cluster age until an orbital radius 
of 13 AU. Since other dust removal processes like 
stellar winds potentially operate at the same time, 
the derived timescale is in fact an upper limit. 

This result thus suggests that at an age of 30- 
40 Myr, M dwarfs can have warm inner debris 
disks. Given that we know already seven M dwarfs 
that harbor extrasolar planets orbiting close to 
their host stars, this finding reinforces the idea 
that planets form around M stars. With M star de- 
bris disks more frequent than previously thought, 
chance occurrences like catastrophic collisions of 
big planetesimals may not have to be invoked for 
their explanation, ruling out the remaining third 
possibility. 

It is interesting to note that a blackbody with 
a peak of its SED at 24 /zm has a temperature of 
120 K, not far from the temperature of w ater sub- 
limat ion, 153 K, defining the "snow line" ( Havashil 



19811 ) , separating the inner region of rocky planet 



formation from the outer region of icy planet for- 
mation and playing a role in the definition of 
"habitable zones" . Ill For M stars, the habit- 
able zone lies at orbital radius of about 0.1 AU 
(pcalo et al.l 120071 ). Indeed, seven M dwarfs are 



now known to harbor one or several extrasolar 
planets at orbital radii o f 0.02 to 2.3 AU (most 
recently Endl et alJ 120081 . see references therein) , 
with masses spanning 0.016-2.1 Mj. The five 
planets with the smallest semi- major axes (< 
0.1 AU) have masses of < 0.07 Mj. During the 
early evolution of an M dwarf, the snow line moves 
inwards due to the decreasing luminosity of the 
central star, enabling the formation of icy proto- 
planets within the first 1 Myr at orbital radii of 
1-4 AU, not taking into acco unt subsequent mi- 
gration (jKennedv et al.l 120071 ). The snow lines of 



1 We can estimate a permissable range of dust temperatures 
for the dust we detect at 24 /im. Comparing the ratio of 
an assumed 3<r excess at 8 )im and the observed 24 /im ex- 
cess with that expected for blackbody emission suggests an 
approximate upper limit of ~ 260 K for a typical source 
without an 8 /an excess. As a strict lower limit, we note 
that the M-dwarf debris disk around AU Mic which does 
not exhibit 24 (im excess em ission has a dust temperature 
of ~ 40 K dRiaz et al.ll2006t) . 



M dwarfs in NGC 2547 have already moved further 
in when compared to their initial location, so that 
rocky protoplanets may have previously formed at 
orbital radii larger than the current snow line. 

Given that some of the few previously known M 
dwarf debris disks have been found at submillime- 
ter wavelengths, it is important to keep in mind 
that material traced in the far infrared is quite 
different from material traced in the submillimc- 
ter radio range. The latter corresponds to mate- 
rial at temperatures of < 2 K at orbital radii o f 
several hundred AU (e.g., iLestrade et al.1 120061 ). 
Thus, different regions are probed by far-infrared 
and submillimeter observations and they cannot 
be directly compared. 

We show the photometric spectral energy dis- 
tributions of the five M dwarfs with the strongest 
excess emission at 24 /im in Fig. [8] Their 24 /xm 
emission is about an order of magnitude stronger 
than the expected photospheric flux. To obtain 
approximate physical properties from the crudely 
constrained photome tric SEDs, we performe d fits 
using disk models bv iRobitaille et al.l ( 20071 ). For 
one of our best-constrained disks, source 23, the 
fits suggest that slightly more than 1% of the to- 
tal luminosity is due to the disk, even though the 
central object is a star with a luminosity of only 
L = 0.028 Lq. To within an order of magnitude, 
the resulting dust disk masses are ~ 10~ 9 A/©, 
corresponding to only a few percent of one lunar 
mass, but in some cases are barely constrained. 
As a major source of uncertainty, these mass es- 
timates scale with the assumed dust properties 
which would be much better constrained with 
(sub-) millimeter data. For all but the densest 
regions of the disk, the model assumes small dust 
particles slightly larger than those in the diffuse 
ISM, and therefore, the mass estimates are a lower 
limit. Reviewing the submillimeter detect ions of 
M dwarf debris disks. ILestrade et al.l (|2006l ) quote 
masses of one to 13 lunar masses for three de- 
bris disks i n the MO range. In marked contrast, 



Chen et al.l (J2005T ) use Spitzer-MIPS data to de- 



rive a minimum dust mass for the debris disk 
of AU Mic of only 10 -4 lunar masses, assuming 
a particle size of 0.2 [im. This mass difference 
may indicate that submillimeter and far-infrared- 
dctcctcd debris disks arc indeed quite different. 
However, a comparison of both methods towards 
the same objects would be needed to confirm this 
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conjecture. 

All of the M dwarf members of NGC 2547 can 
be assumed to be of the cluster age, 30-40 Myr, 
possibly indicating that at this age, the evolu- 
tionary stage of M dwarfs makes it easier to de- 
tect such debris disks. For debris disks around 
stars of earlier spectral types, there have already 
been several studies concerning age trends. In a 
study of ma i n-sequ ence A stars older than 20 Myr, 
Rieke et al.l ( 20051 ) find that the magnitude of ex- 
cess at 24 /mi declines with a to/t dependence, 
with to w 150 Myr. Studying the double clus- 
ter h and x Persei i n context with this result and 
other observations, ICurrie et alJ ( 20081 ) conclude 
that the debris disk excess emission of A stars 
peaks at about the age of this cluster, 13 Myr. 
On lo nger timescale s, the evolution is less clear. 
While ITrilling et ail (|2008l ) did not find a clear 
age trend studying FIR excess emission of solar- 
type fiel d stars with ages ran ging from 200 Myr to 
10 Gyr, iMever et alJ (|2008j ). also studying solar- 
type stars, found an excess rate decreasing from 
18% at ages 3-30 Myr t o 2% at ages 0.3-3 Gyr. 
Hillenbrand et all (J2008I ) find that the 70 /mi ex- 
cess emission of 328 FGK stars, tracing cooler dust 
compared to 24 /mi, peaks at ages of 30-100 Myr. 

For a solar-mass star, the debris for rocky 
planet formation peaks after a few million years 
for m aterial at 1 AU (JKenvon fe Bromlevl I2004L 
2005) . To scale this result for an M dwarf, we can 
use a simple relation for the disk mass, A/disk oc 
M Rt , a , ri and the collision ti mescale, i co u tx M t 



-1/2 



(JKenvon fc Bromlevl 120081 ). So, if the M dwarf is 
three times less massive than the Sun, we would 
expect planet formation to take about five times 
longer around an M dwarf, i.e., 10-20 Myr, not 
taking into account minor complications due to 
the moving snow line. This indicates that the M 
dwarfs in NGC 2547 are beyond but still close to 
the maximum predicted debris production. Espe- 
cially with regard to stars of earlier spectral types, 
it is important to keep in mind that this argument 
only allows to compare the evolutionary timescale 
of debris at an orbital radius of 1 AU. For A stars, 
the radius probed by 24 /iin observations is beyond 
20 AU, and the timescales at these radii may be 
very different. 

Currently, the early debris-disk evolution of M 
dwarfs is poorly constrained from an obse rvational 
point of view. IWeinberger et al.l ( 2004 ) did not 



find significant mid-infrared excess (12 and 18 /mi) 
towards 16 members of the young, nearby TW 
Hydrae association and conclude that apparently 
any planet formation in the terrestrial planet re- 
gion was rapidly completed. In a study target- 
ing stars at age simila r to the age of NGC 2547, 
Mamaiek et alJ ( 20041 ) observed members of the 



30-Myr-old Tucana-Horologium Association, in- 
cluding two M dwarfs, at a wavelength of ~10 /im 
(JV band). They do not find excess emission to- 
wards the two M dwarfs, and their quoted uncer- 
tainties would allow excess emission at only ~10% 
of the expected photospheric flux level. Look- 
ing at our SEDs in Fig. [51 if the excess emis- 
sion is modeled correctly, it would be far more 
difficul t to detect at 1 /mi. The M dwarfs stud- 
ied by IPlavchan et al.l (|2005l ). and not found to 
exhibit excess emission at 11.7 /im, include one 
young source and nine more objects with ages of 
>600 Myr. The single young star is GJ 3305, 
a member of t he Pic movi n g grou p at an age 



of -12 Myr. iGautier et all (|2007l ) found that 



field M dwarfs older than 1 Gyr do not show far- 
infrared excess emission. Due to the nearby loca- 
tion of these sources (d < 20 pc), the absolute sen- 
sitivity of these observations is much higher than 
that of our dataset. This upper limit for the evolu- 
tionary timescale of M dwarf debris disks indicates 
that planet formation around M dwarfs has ceased 
by this time. 

We finally consider the excess fr actions of stars 
with e arlier spectral types. While iGorlova et al. 
(|2007l ). restricting the discussion to their com- 
pleteness limit, found 30 - 45% of the B-F mem- 
bers to show excess emission at 24 /mi, we inter- 
estingly find virtually no G star excess emission 
at this wavelength. As pointed out above, while 
we find excess emission towards K stars, we may 
not be complete throughout their R-K range, and 
many of the "highly probable" K members may be 
background giants. Even so, the 24 /mi excess rate 
for both G and K stars, using our conservative cri- 
terion for excess sources, is at most a few percent. 
In th e most compr ehensive sample of field stars to 
date, ITrilling et al.1 (|2008) find statistically indis- 
tinguishable excess rates for A, F, G, and K stars, 
with an average age of the solar-type sample of 
5 Gyr. The fact that we do not find excess emis- 
sion towards G stars can possibly be understood 
using an argument similar to the one used in our 
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discussion of M dwarf debris disks: given the clus- 
ter age, the solar-type stars are far beyond their 
peak in dust production. In marked contrast, the 
excess rate for M dwarfs is surprisingly high, even 
though it is a lower limit for two different reasons: 
first, our 24 /im detections of M dwarfs are not 
complete and second, we use a conservative crite- 
rion to define excess sources. It remains unclear, 
however, to what degree the excess rates of G, K, 
and M stars are skewed due to contamination by 
non-members. 

6. Summary 

In an exceptionally deep 24 /im observation 
of NGC 2547, we detect eleven M dwarf excess 
sources interpreted as debris disks, nine of which 
are new. These constitute more than 42% of 
all significant excess sources that we find across 
the range from B to M stars. The eight M 
dwarf excess sources that were covered by IRAC 
observations do not show significant excess in 
the se mid-infrared ban ds (except for source 23, 
see lGorlova et alJl2007l ). The dust masses of these 



disks are crudely estimated to be on the order of 
a few percent of a lunar mass. The identification 
as debris disks hinges on the observation that the 
dust removal timescale due to processes like PR 
drag and stellar winds is much smaller than the 
age of the cluster. Thus, there must be a mech- 
anism to replenish the dust particle population, 
e.g., by collisions. Given the cluster age, the ob- 
served M dwarfs appear to be close to, maybe 
beyond their maximum predicted debris produc- 
tion. Interestingly, the observed debris material is 
located close to the snow line and the habitable 
zone, at less than 1 AU from the central stars. 
Since even our deep dataset is not complete for 
M dwarfs, we cannot give a debris disk fraction in 
terms of the total population, but it appears to be 
higher than for G and K stars at the same age. 
In fact, we find a minimum in excess emission in 
the G stars, suggesting that their maximum debris 
production occurs much earlier than at 30-40 Myr. 
With these new discoveries, the number of 
known M star debris disks more than doubles by 
rising from six to fifteen. The previously perceived 
"dearth" of M dwarf debris disks apparently was 
due to small sample sizes both in numbers of ob- 
jects studied and in different ages probed. To bet- 



ter constrain the debris-disk evolution and planet 
formation timescales of M dwarfs, larger samples 
of objects at ages both younger and older than 
NGC 2547 are needed. 
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